A current driven instability in parallel, relativistic 
shocks 



B. Reville^, J.G. Kirk 2 , P. Duffy 1 

1 School of Mathematical Sciences, University College Dublin, Ireland 
2 Max-Planck-Institut fur Kernphysik, 69029 Heidelberg, Germany 

E-mail: Tbrian.reville@ucd.ie 

Abstract. Recently, Bell [2] has reanalysed the problem of wave excitation by 
cosmic rays propagating in the pre-cursor region of a supernova remnant shock 
front. He pointed out a strong, non-resonant, current-driven instability that had been 
overlooked in the kinetic treatments by Achterberg [T] and McKenzie & Volk ^0]> and 
suggested that it is responsible for substantial amplification of the ambient magnetic 
field. Magnetic field amplification is also an important issue in the problem of the 
formation and structure of relativistic shock fronts, particularly in relation to models 
of gamma-ray bursts. We have therefore generalised the linear analysis to apply to this 
case, assuming a relativistic background plasma and a monoenergetic, unidirectional 
incoming proton beam. We find essentially the same non-resonant instability noticed 
by Bell, and show that also under GRB conditions, it grows much faster than the 
resonant waves. We quantify the extent to which thermal effects in the background 
plasma limit the maximum growth rate. 
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1. Introduction 

The acceleration of cosmic rays at the shock front that bounds a supernova remnant 
is thought to proceed via the diffusive first-order Fermi mechanism (for reviews, see 
jniEHE]) In the shock precursor, where the cosmic rays stream at roughly the shock 
speed, Alfven waves can grow as a result of interacting at the cosmic ray cyclotron 
resonance. In turn, these waves provide the pitch-angle scattering essential for the 
acceleration process Using a kinetic approach, Achterberg l.j found that the 

presence of cosmic rays has no significant impact on the nature of the plasma modes. 
However, Bell using a hybrid MHD/kinetic theory approach, has recently challenged 
this conclusion. He identifies a low frequency, non- resonant, electromagnetic mode that 
can be strongly driven by the current induced in the plasma by the streaming cosmic 
rays. Simple analytic estimates and MHD simulations both suggest that the nonlinear 
evolution of this instability leads to very strong amplification of the ambient magnetic 
field. In turn, this may result in a cosmic-ray acceleration rate that is much more rapid 
than previously thought [3J Ej . 

Gamma-ray bursts drive highly relativistic outflows with Lorentz factors of 100 
or more ^3]. On interacting with the surrounding medium, a shock front forms, 
but the mechanism by which this happens is controversial [18j . Observations of 
the afterglows of gamma-ray bursts suggest that the ambient magnetic field must 
be amplified substantially at the shock front. The relativistic Weibel or two-stream 
instability, that generates a small-scale magnetic field in a previously unmagnetized 
plasma, has been investigated in this connection [TT], but appears to saturate at a 
relatively low amplitude [3 GH|. In recent years, particle in cell (P.I.C.) simulations 
have been used in studies of collisionless shocks [IB], in particular with reference to the 
Weibel instability [7] US] Ej. The instability that is presented in this paper has not 
been identified to date. However, relativistic P.I.C. electron-proton simulations in three 
dimensions still present a major challenge, even to the best computational resources 
currently available. 

When two electron-proton plasmas moving at high speed relative to each other 
collide and interpenetrate, the first process that takes place is the equilibration of 
the electrons [Ej. However, this releases only a small fraction of the available free 
energy. The shock front responsible for the thermalisation of the bulk of the energy is 
mediated by the interaction of the two counterstreaming proton plasmas. Thus, Pohl 
& Schlickeiser ^3] attacked the shock formation problem starting with a monoenergetic 
relativistic proton beam penetrating a relatively dense, cool background in a direction 
aligned with the ambient magnetic field. They then computed the isotropisation rate 
that results from the resonant excitation of Alfven waves in the background plasma. 
Here, we use the same initial conditions, but, because this physical situation resembles 
that in the precursor region of a SNR shock front, we analyse instead the growth rate 
of the relativistic analogue of the non-resonant mode discovered by Bell [2 . We find it 
grows much faster than the resonant mode, and, as in the case of SNR, can be expected 



A current driven instability in relativistic shocks 



3 



to generate a substantial magnetic field transverse to the beam direction. However, we 
find the instability is quite sensitive to damping by thermal effects, once the background 
plasma is heated. 

2. Dispersion relation 

The linear dispersion relation for the propagation of transverse waves parallel to the 
magnetic field in a plasma made up of components labelled by j is: 

f i-^(M) = o (l) 

j 

where m = ck/u is the refractive index and Xj{k,uj) is the susceptibility of the j'th 
component, for wavenumber k and frequency u (> 0). 
The three components in the case we consider are: 

(i) Protons that make up the background or downstream distribution, denoted by 
j = p. These have a number density n p and a thermal distribution with 
(dimensionless) temperature 9 P = k B T p /m p c 2 . The waves will be analysed in 
the lab. frame, chosen such that these particles have a vanishing net drift speed 
c/3 p = 0. 

(ii) Electrons that stem partly from the downstream plasma and partly from the 
incoming upstream plasma or beam. They also have a thermal distribution with 
lab. frame density n e and temperature e . Their drift speed in the lab. frame is 
c(3 e and their Lorentz factor r e = (1 — 0£) 1 ^ 2 . 

(iii) Protons of the upstream medium that form a monoenergetic, unidirectional 
incoming beam along the magnetic field. Their lab. frame density is and their 
drift Lorenz factor T^. The beam distribution function is 

f b (u) = (— ^) 5(u±)5(u l{ - T b (3 b ) (2) 

Following Achterberg Ij, we impose the conditions of overall charge neutrality and 
zero net current on these components. This implies: 

I>PiM? =° ( 3 ) 
j 

and 

3 

where u P j = (Aitrijqlj/mj} ^ , denotes the plasma frequency of the j'th component that 
consists of particles of charge q 3 - and mass rrij, and u C j = qjB/irijC, denotes its cyclotron 
frequency. In the present case m b = m p is the proton mass and qb = % = —q e = e is 
the modulus of the electronic charge. 
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For parallel propagation, the susceptibility Xj is given by Yoon (201 • The quantity 
uj 2 Xj, which is invariant to Lorentz boosts along the magnetic field, can be written: 

-ivy + cku\\ u\ (c 2 k 2 — u 2 



d 3 u 



7 



D («||) 2 D 2 



(4) 



Here, fj(u) is the distribution function of particles of four velocity c(j,ju), normalised 



such that / d ufj = 1, the components of u parallel and perpendicular to the magnetic 
field direction are u\\ and u±, and the resonant denominator is 

D(uA = suj C3 (l + Z (uA) (5) 



(6) 



with 

uj'-f — cku\\ 



eu} cj 



and 7 = yT + u\ + u 2 . The waves are circularly polarised, with e — 1 corresponding to 
left-handed, and e — — 1 to right-handed waves, for k > 0. 

The waves we consider are non-resonant for the electrons and background protons. 
Furthermore, these components are "magnetized", in the sense that for all relevant 
values of u\\, the resonant denominator defined in Eq. (jHJ) can be expanded using as a 
small parameter Z (u\\J <C 1: 



D n (u\\) lew. 



I I I . „ ! \ I' \ II I / „.., / \ \ 

Inserting this expansion into Eq. (jlj), one finds, using / d 3 ufj {u)u\\/^ = (3j, for the 
zeroth x an d first order x contributions: 

u^xf (k, u,) = (cfrk -u)-^L (Jk 2 - u 2 ) Ulh) . (8) 



u \f ( k , u)=J (u 2 - look («||) . + c 2 k 2 (u\h). 

(c 2 k 2 - cu 2 ) (ck (ulu\\h). - co (u 2 ± )^j (9) 



where (. . .) j = J d 3 u . . . fj (u). 

The susceptibilities Xj describe the currents induced by the wave field in each 
plasma component. The first term on the RHS of Eq. (JHJ) is proportional to the charge 
and current density of the j'th component in the unperturbed state. In a neutral, 
current-free plasma in which all components are magnetized, i.e., can be treated using 
the expansion (jZJ), these terms cancel, according to Eqs (jHJ) and (J1J). However, in the 
plasma we consider, only the background protons and electrons are magnetized, and 
these terms do not cancel: 



ujL , _ . . ^ cu 2 



E ^xf (K w ) = -J*> - cp b k) - £ -4 U 2 k 2 - u 2 ) (ulh) . (10) 



j=p,e £UJ cb j=p,e "~cj 




Figure 1. Dispersion curves for the illustrative parameter set: va = 2 x 10~ 5 c, 
n>b = "-p/3, Tb = 10 for two temperatures: P = 10 _1 (right panel) and 10 -3 (left 
panel). The real part of the frequency of each mode is shown as a dashed line, the 
imaginary part as a solid line. The left(right)-handed waves are shown in grey(black). 



We describe the waves in a reference frame in which the background protons have 
zero drift, and assume they have an isotropic distribution in this frame. In this case 
(«ll) p = (u 2 ± u ]l /j) p = and (u\h) p = 2 («J/t) p - 

Inserting the distribution function (j2J) into Eq.(J3J), we find the susceptibility of the 
beam protons is 

LU X h{k,0o) = . r (11) 

EU) ch - T b (c(J h k - LO) 

If TbC/3b \k\ <C |cu c b|, the expansion ((Zj) can also be used for the beam protons, so that 
the plasma is fully compensated. However, modes of shorter wavelength, such that 

|^cp,e| /@p,c > C\k\ > |w cb |/r b (12) 

are unable to induce a compensating current in the beam particles. The overall 
susceptibility can then be written: 

^X~— — — -, + —5- + [c k - w ) (u ± ) 13 

eto c eu c + uj' v A ecu* v / \ / p 

where we have neglected the electron response, except for its contribution to the overall 
current in the first term and to the quantity v a, defined as the corrected non-relativistic 
expression for the speed of an Alfven wave in the background plasma |2()j : 



2 2 
V A = C 



UJ 2 



. p,e ^cj 



E 7X <7>, + </ 2 7 



-i 



(14) 



The quantity u' = r b (uj — /3 b cfc) denotes the wave frequency as seen in the frame 
of the beam particles, and the plasma frequency of the beam particles in this frame is 
given by <j\ = ul/T b . 
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3. Wave modes 



For a cold background plasma, and for low frequency modes such that u' « — T^Pbck 
with |u/| ^> |^cb|, one recovers from Eq. (fT^j) a simple form analogous to that discussed 
by Bell 0: 

u 2 = (c 2 k 2 + (15) 

which gives a purely growing mode (Re (a;) = 0) for the right-handed polarisation 
e = — 1, provided the driving by the beam- induced current is sufficiently strong. For 
v a <C c, the mode reaches a maximum growth rate 

Im(ct>) = -— PbUpp (16) 

which is independent of the magnetic field strength. The corresponding wave number is 

*W = (17) 
2n p v A 

The (non-resonant) thermal effects on this mode arising from the term containing 
(u 2 ± ) are easily analysed in the case of weak magnetic field v& <C c, since then u 2 <C c 2 k 2 
and u/ ~ — Tbf3\ } ck. The dispersion relation is then approximately 

c 2 \ ew 2 £T b /3 b < / 

where dimensionless units are used in which k = T\)f3\ } ck/uj c and Co = T^P^u / 'u c . The 
e = —1 mode reaches a maximum growth rate of 

In terms of the dimensionless temperature p of the background protons for the 
appropriately normalised Jiittner-Synge distribution 



u 



I) = 166J [1 + 3(1 + P )] e- 1 ^/ [K 2 (1/6 P )] (20) 

so that the growth rate given by Eq. ()19j) falls off rapidly with increasing temperature. 

After multiplication by the factor eu c — ou', Eq. (0) together with Eq. ()13|) gives 
the dispersion relation as a fourth order polynomial in u with real coefficients. Keeping 
all terms in this expression, which is valid for wavenumbers k <C T^P^/Q, far from 
resonance with the background protons, the low frequency modes are shown in Fig. ^ 



4. Discussion 

From Fig. ^ it can be seen that the e = 1 mode (grey curves) resonates with the beam 
particles when k is close to their inverse gyro-radius, i.e., k ~ 1. This is the cyclotron 
resonant mode considered by, amongst others, Pohl & Schlickeiser [15j. However, for 
low to moderate values of the background temperature, its growth rate is very much 
slower than that of the non-resonant, current-driven mode (black curves). Near the 
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peak of the unstable band of the non-resonant mode, the waves are strongly modified 
by the beam and do not resemble Alfven waves: the phase velocity of the non-resonant 
mode is lower than that of the resonant mode, and both speeds generally exceed va 
substantially. At very large wavelengths, the modes converge and become Alfven waves. 
This also happens for large k values outside the unstable band, but only if these do not 
resonate with the background protons. 

Eqs ()16p and (|17[) . show that in a cold background with ~ n p , the non-resonant 
mode has a maximum growth rate comparable to that of electrostatic oscillations, 
although at much shorter (parallel) wavelengths. Only a relatively modest temperature 
6 P ~ cg> c /cg> pp suffices to damp such waves, as can be seen from Eq. (fTH|) . Physically, this 
can be understood by noting that, as the background temperature rises, an increasing 
number of particles react to the wave perturbation as if they were unmagnetized. This 
reduces the effective current driving the instability. Ultimately, if all particles react in 
the same way to the perturbation (i.e., all are magnetized or all are unmagnetized) the 
plasma behaves as a fully compensated system with no driving current. Eq. (fTH|) also 
shows that, for fixed temperature, the growth rate scales with the cyclotron frequency 
of the background protons. For P « 1 and ~ n p , the maximum growth rate is 
approximately equal to uj c . 

The crucial question of the nonlinear evolution of the system has been discussed in 
the nonrelativistic case by Bell [2] and Milosavljevic & Nakar [12] ■ Saturation can be 
expected when the currents associated with the growing waves become comparable to 
the current induced by the beam. Applying this to the relativistic case, one finds for 
the amplitude B w of the magnetic field of the perturbations 

k A B w w q h n h cp (21) 

If we assume that the largest wavelength in our system is of the order of the gyroradius 
of the beam protons, we can estimate the maximum magnetic field energy density to be 

S « l -n h Tm v c 2 (22) 

This suggests that the entire energy content of the beam is, at least initially, 
transferred into transverse field fluctuations. The efficiency of magnetic field generation 
by this mechanism (i.e., the ratio of the saturated field energy density to that in the 
incoming stream) is somewhat higher than that attributed to the Weibel instability 
in the case of a relativistic shock in a pair plasma ^H]- The Weibel instability is not 
thought to be effective in mediating shocks in electron/proton plasmas [5]. However, the 
efficacy of the instability described here depends strongly on the ability of the plasma 
to convert the energy input at small scales into a larger scale magnetic field. 
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